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Abstract 



In order to determine (as best as possible) the nature of the dark matter particle (mass and decoupling temperature) we compare 
jhe theoretical evolution of density fluctuations computed from first principles since the end of inflation till today to the observed 
properties of galaxies as the effective core density, core radius and surface density. We match the theoretically computed surface 
Q^idensity to its observed value in order to obtain: (i) the decreasing of the phase-space density since equilibration till today (ii) the 
mass of the dark matter particle and the decoupling temperature (iii) the type of density profile (core or cusp). The dark matter 
^ particle mass turns to be between 1 and 2 keV and the decoupling temperature turns to be above 100 GeV. keV dark matter 
particles necessarily produce cored density profiles while wimps (m ~ 100 GeV, ~ 5 GeV) inevitably produce cusped profiles 



' at scales about 0.003 pc. We compute in addition the halo radius ro, the halo central density po and the halo velocity -sf^haio- they 
all reproduce the observed values within one order of magnitude. These results are independent of the particle model and vary 
Very little with the statistics of the dark matter particle. Moreover, our results for the surface density using typical CDM wimps 
agree with the values obtained from CDM simulations (which are five orders of magnitude larger than the surface density observed 
O-iValues). DM particles with mass in the keV scale well reproduce the observed properties of galaxies as the effective core density, 
Q cored profiles and the value of the surface density. The linear framework presented here turns to apply both to keV and to GeV 
"mass scale DM particles. keV scale particles reproduce all observed galaxy magnitudes within one order of magnitude while GeV 
mass particles disagree with observations in up to eleven orders of magnitude. 
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1. Introduction 



Since several years and more recently H [il H B it 



X 



has been stressed that basic galaxy parameters as mass, size, 
baryon-fraction, central density, are not independent from each 
pther but in fact all of them do depend on one parameter that 
works as a galaxy identifier In fact there exist functional re- 
lations that constrain the different galaxy parameters in such a 
way that the galaxy structure depends essentially on one param- 
eter ( Bill and references therein). 

' These functional relations may play for galaxies the role that 
J:he equations of state play in thermodynamical systems. 

First, let us remind that the density of DM halos around 
galaxies is usually well reproduced by means of dark halos with 
a cored distribution (author?) Isllsjl, where ro is the core ra- 
dius, po is the central density hm p(r) - po and p(r) for r < ro 

r— >0 

is approximately constant. Recent findings highlight the quan- 
tity yUo = rQ Po proportional to the halo central surface density 
defined as 



Jo 



p(0, 0, X3) iix3 where r = (xi, X2, jcj) 



galactic systems (spirals, dwarf irregular and spheroidals, ellip- 
tics) spanning over 14 magnitudes in luminosity and over dif- 
ferent Hubble types. More precisely, all galaxies seem to have 
the same value for //o, namely jiQ — 120 Mo/pc 
It is remarkable that at the same time other important structural 
quantities as ro, po, the baryon-fraction and the galaxy mass 
vary orders of magnitude from one galaxy to another 

The constancy of po is unlikely to be a coincidence and prob- 
ably has a deep physical meaning in the dynamics of galaxy 
formation. It must be stressed that po is the only dimensionful 
quantity which is constant among galaxies. 

By analogy with the theory of phase transitions in statisti- 
cal physics we find useful to call 'universal' those quantities 
which take the same value for a large set of galaxies while non- 
universal quantities vary orders of magnitude from one galaxy 
to another. In this context the quantities called universal take the 
same value up to ±20% for different galaxies. That is, they are 
only approximately 'universal'. Other known universal quan- 
tity in the above sense is the shape of the density profile when 
expressed as a function of r/ro and normalized to unit at r = 0. 

We use here the linear approximation to the Boltzmann- 



Vlasov equation 11131 12511 to compute the DM density profile 



where X3 goes along the line of sight. The quantity po is found 
nearly constant and independent of luminosity in different 



and from it the surface density po. 

In this paper, from first principles and using the primordial 
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power spectrum we follow the evolution of the gravitational 
collapse of a perturbation of mass M ~ 3 x lO'^M© and de- 
rive the resulting linear halo density profile. This reproduces 
the phase of fast accretion found in A^-body simulations. As a 
result we obtain robust predictions for the properties of DM ha- 
los from first principles. In the case of ACDM our results agree 
with the A^-body ACDM simulations and in the case of AWDM 
our results agree with the observations. 

We obtain a very good fit of the linear profile to the Burkert 
profile which determines the relation between ro and the free- 
streaming length. This shows that our linear profile can appro- 
priately fit halo density data. The linear approximation to the 
Boltzmann-Vlasov equation improves for large halo radius ro 
(see fig. U}- Since universal quantities take the same value for 
large or small galaxies one can expect good results from the lin- 
ear approximation which improves for large and less compact 
galaxies. 

We also compute non-universal galaxy quantities as the halo 
radius, galaxy mass, halo central density and squared halo ve- 
locity. We find that the linear approximation provides halo cen- 
tral densities smaller than or in the range of the observations, 
and halo velocities larger than the observed ones by a factor 
between 1 and 10. Notice that the determination of these non- 
universal galaxy quantities is unrelated to the computation of 
the DM particle mass. Hence, our linear calculations amply 
suffice for our specific goals to determine properties of the DM 
particles. 

Clearly, for non-universal quantities one has to include the 
nonlinear effects (including for instance mergers) to exactly re- 
produce the observations. Notice that general arguments based 
on the Boltzmann-Vlasov equation show that cores cannot be- 
come cusps through mixing and mergers and that cusps do not 
become steeper neither shallower through mixing and mergers 

Moreover, recent A^-body ACDM simulations (Acquarius) 
have found that the DM halos form in a sort of "monolithical" 
way. Their inner regions, that contain the visible galaxies, are 
found to be stable since early times and contrary to previous 
believes, major mergers (i.e. those with progenitor mass ra- 
tios greater than 1 : 10) are found to contribute little to their total 
mass growth 1391]. This indicates that nonlinearities (i.e. merg- 
ers) have a reduced importance. Minor mergers, secondary in- 
fall, rare major mergers are certainly important for details, but 
the gross features of DM halos are determined during the fast- 
accretion phase of their gravitational collapse, as the history of 
the quasar-galaxies coevolution also seems to indicate ll20ll . 

The halo formation essentially consists of two main phases: 
A first fast accretion phase (that can be treated by the linear 
approximation), and a second subsequent slow accretion phase 
with mergers and infalls, that have a random character and that 
can only be described by numerical simulations. This second 
phase does not have an essential influence in the shape of the 
halo profile. Thus, in order to explain the observed halo profiles 
one just needs to describe the first phase of halo formation, as 
we do here in this paper. 



Therefore, the cored and cusped character we find for the 
linear profiles depending on the DM particle mass considered 
(keV and GeV mass scale, respectively) should remain valid 
after mixing and mergers are taken into account. 

Clearly, to build the theory of galaxy formation simulations 
are needed using the correct linear primordial power spectrum 
as input. Our work is only indicative in that respect. 

We combine in this paper observed properties of galaxies as 
the effective core density and the core radius with the theoreti- 
cal linear evolution of density fluctuations computed from first 
principles since the end of inflation till today. We compute the 
galaxy surface density and match it with the observed values. 
Within the same scheme, we derive analytically the halo radius 
ro and the factor Z characterizing the reduction of the phase- 
space density since equilibration till today. In this way, we find 
the value of the mass of the dark matter particle which turns 
to be between 1 and 2 keV, and the number of ultrarelativistic 
degrees of freedom of the dark matter coupling at decoupling 
gd, or similarly, the decoupling temperature Tj which turns to 
be above 100 GeV. Evidence based on the phase space density 
pointing towards a DM particle mass in the keV scale was pre- 
sented in refs. ilsi lToli . 

We consider in this paper the whole range of galaxy virial 
masses going from 5 to 300 xlO" Mq. 

Non-universal quantities are reproduced in the linear approx- 
imation within a factor of ten. 

The analytic expressions we derived for the density pro- 
file, and the mass of the dark matter particle also imply that 
keV dark matter particles always produce cored density profiles 
while heavy dark matter particles as wimps {m =100 GeV, Td 
= 5 GeV) inevitably produce cusped profiles at scales of 0.03 
pc. These results are independent of the particle model, vary 
very little with the statistics of the dark matter particle and are 
all obtained in the linear framework. 

Analytic methods have been used since some time to de- 
rive galaxy properties using the primordial power of the density 
fluctuations [see for example |21, 2^] and using the spherical 
model (author?) H 



The theoretical treatment presented here captures many es- 
sential features of dark matter, allowing to determine its nature. 
Our treatment also applies to CDM: if we use the CDM surface 
density value obtained from CDM simulations ||22|. we deter- 
mine [sec. |8| a dark matter particle mass in the wimps mass 
scale (GeV), fully consistent with CDM simulations. 

This paper is organized as follows: Sec. 2 presents galaxy 
data and empirical formulas relating basic galaxy parameters; 
sec. 3 deals with the phase-space density; sec. 4 contains 
our theoretical results for the density profile from the linearized 
Boltzmann-Vlasov equation. In sec. 5 we derive the DM par- 
ticle mass and the decoupling temperature from the theoretical 
and observed galaxy surface density, in sec. 6 we compute non- 
universal galaxy properties and in sec. 7 we derive the profiles 
for keV scale DM particles and for wimps (cored vs. cusped 
profiles). In sec. 8 we present our conclusions. 
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the surface density is a constant over a large number of galaxies 
of different kinds. 

Notice that the surface density of ordinary matter in lumi- 
nous galaxies is about a factor 4 larger than the surface density 
value for dark matter UtIi . Clusters of galaxies, exhibit a dark 
matter surface density about a factor 4 or 5 times larger than 
that of dwarf, elliptical and spiral galaxies JsKIstl- Such differ- 
ence could be due to a baryons effect, which study is beyond 
the scope of this paper For clusters of galaxies ro is 4 to 50 
times larger than for the galaxies in Table 1 and the masses are 
100 to 4000 larger than the masses of the galaxies in Table 1. 
Namely, the variation of //o from galaxies to clusters of galax- 
ies is a much small factor than the change in ro and in the total 
mass. We choose for the present work the data from galaxies 
in Table 1 (further discussion on clusters of galaxies is given in 
sec. 6). 



Figure 1: The common logarithm of the observed surface density fJoobs in 
(MeV)^/{rf vs. the common logarithm of the core radius ro in kpc. Notice 
that in galaxies both rg and po vary by a factor of thousand while /iq varies only 
by about ±20%. 



2. DM halos around galaxies: the observational scenario 

The kinematics of about several thousands disk galaxies, 
described by the Universal Rotation Curves of Spirals, and 
the information obtained from other tracers of the gravita- 
tional field of galaxies, including the dispersion velocities of 
spheroidals and the weak lensing measurements ( 113 ill and ref- 
erences therein) found that the density of the dark matter ha- 
los around galaxies of different kinds, different luminosity and 
Hubble types is well represented, out to the galaxy virial radius, 
by a Burkert profile 



p(r) = Po Fb\—] , Fb(x) = ^- , x = — 

FO B^^J (l+x)(l+x2) ro ^ 

where po stands for the effective core density and ro for the core 
radius. The Burkert profile satisfactorily fits the astronomical 
observations and we use the observed values of po vs. ro for 
DM dominated spiral galaxies given in llsill . 

The structural halo parameters po and ro are found to be re- 
lated, it is worth to compute from them the virial mass M,.,-,- in 
terms of the core radius ro ( 113 111 and references therein) 



M,i,. 



■iai 



0.320 ( — 
\kpc 



101 'Mo 

The surface density jUo is defined as: 

= Po ro 



(2) 



(3) 



We display in Table 1 the values of the observed surface den- 
sity /uoobs in (MeV)-'/(fi,^ c'') and the corresponding core ra- 
dius ro. We plot in fig. [T]the observed surface density fiQobs in 
(MeV)-'/(?i^ c^) vs. the core radius ro. 

Notice that in galaxies both ro and po vary by a factor lO-' 
while yUo varies only by less than ±20%. 5 kpc < rp < 50 kpc 



for normal spiral galaxies. Therefore, as stressed by 11261114134 1 



3. The invariant phase-space density of DM galaxy halos 

The invariant phase-space density is defined by |23, Tllslfiol] 
G = 4 (4) 
where 

1 

(5) 



= - <v > 

is the velocity dispersion. Q is invariant under the expansion of 
the universe and decreases due to self-gravity interactions Issll 
from its primordial value Qp to the volume average value Qhaio 
of the galaxies today: 



Qbalo — 2 



where 



halo — 



Phalo 



Ppn 



(6) 



(7) 



halo 



This equation defines the factor Z ifioll . Z is larger than unity 
and its value depends on the galaxy considered. 

Notice that Qp only depends on the properties of the DM 
particle and its primordial distribution function [see eq.(fT2li be- 
low]. 

The squared halo velocity v^^,^(r) follows from the virial the- 
orem combined with the Burkert profile eq.O with the result 
ill 



^hah 



(r) = 2nG 



porl 



1 2 

ln(l + x) - arctan.;c + - ln(l + x ) 



r 
ro 



Qhaio follows by averaging p(r) and v^^j^(r) over the volume 
using the density itself p(r) as weight function (see Appendix 
Appendix Ai. From eqs. (HJ, ©, and (O we obtain [see 



eq.( IA.4l )l. 

Qhaio - 



0.069 

G^- yfPor^ 



(9) 



/"n in knc 


/vn in MeV"^ 






0. i 


U.04 iU 


7 Q 




1 n 7 
i w.z 


OA? 1 n4 


13.3 


0.61 10^ 


17.3 


0.60 10"* 


22.6 


0.60 lO'* 


29.4 


0.59 lO'* 


38.7 


0.57 10-* 


51.8 


0.55 10-* 



Table 1: The observed core radius ro and the observed surface density fioobs ■ 
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Figure 2: The common logaiithm of the phase-space density Qi,aio obtained 
from eq.^9) using the data in Table 1 vs. the virial mass of the galaxy niy = 
A^v;rra//[10"Mo]. [niy is related to the data ro in Table 1 through eq.^2)]- 



For a NFW profile. 



p(r) = 



r / r 
- 1 + - 



2 ' 



(10) 



we get for Qhaio following the same steps as for the Burkert 
profile 



halo 



0.324 



(11) 



Both results eqs. ^ and (fTTT i are of the same order of mag- 
nitude and differ by a factor ~ 5. Since Q ~ m'^ as shown 
below in eq. (fT2l i. using the cuspy NFW profile instead of the 
cored Burkert profile only may change the DM particle mass 
by a factor ~ 1 .5 keeping its order of magnitude. 



We plot in fig. |2]the common logarithm of the phase-space 
density Qi,aio vs. the virial mass of the galaxy ot,>. Notice that 
the virial mass of the galaxy m,, is related to the halo radius ro 
through eq.©. 

The primordial invariant phase-space density Qj, can be eval- 
uated in the radiation dominated (RD) era with the result ifioll 



3 V3 
2l^ 



m 



(12) 



where h and h are the dimensionless momenta of the particle 
DM primordial distribution function, g is the number of inter- 
nal degrees of freedom of the DM particle {g - 2 for Dirac 
fermions). For example, for Dirac fermions of mass m that de- 
coupled ultrarelativistically at thermal equilibrium we have. 



Qp = 0.020395 . 



(13) 



Similar expressions and values are obtained for bosons and for 
particles decoupling ultrarelativistically out of thermal equilib- 
rium [10]. 

The covariant decoupling temperature can be expressed 
in terms of the number of ultrarelativistic degrees of freedom at 
decoupling gd by using entropy conservation iQ] 



(14) 



gd can be expressed as ifioll 
2l 



gd 



3i ni ^DM Pi 



i 7^ i 



(15) 



where Ty is the CMB temperature today, Q.dm the DM cosmo- 
logical fraction and pc the critical density of the universe: 



Ty = 0.2348 meV , Q.dm = 0.228 , 



Pc = (2.518 meV)7(?i c") 
here 1 meV = 10"^^ eV. 



(16) 
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4. Theoretical results. The linear Boltzmann-Vlasov equa- 
tion. 

We will now evolve the density fluctuations from the end of 
inflation till today in the standard model of the Universe. This 
evolution provides the phase-space density Qhau, and the sur- 
face density //q today. The density fluctuations follow from the 
distribution function which evolves according to the non-linear 
Boltzmann-Vlasov equation. The evolution is practically linear 
in the RD era and in the MD era before structure formation. 
That is, we can use the linear Boltzmann-Vlasov for redshift 
z > 30. For z < 30 non-linearities are relevant and one should 
use the non-linear Boltzmann-Vlasov equation or, alternatively, 
perform A^-body simulations. 

It must be noticed that the resolution of the linearized 
Boltzmann-Vlasov equation from the end of inflation till today 
provides a good approximated picture of the structures today 
1 1 ill . From this linear evolution of the dark matter fluctuations 



A(^, z) we obtained the linearized density profile pii„(r) 111 1 1. 
We followed the density fluctuations in the RD era according 
to the results in 11311 and 12411 . It is convenient to recast the 
linearized Boltzmann-Vlasov equation in the matter dominated 
(MD) era as an integral equation, the Gilbert equation (l^ . We 
solved the Gilbert equation | H , 3l to obtain the density fluctu- 
ations A{k, z) till today with the result 



Bose-Einstein 
Fermi-Dirac 
Maxwell-Boltzmann 




A{k,z)'^ - T(k){l+z,n)A{k,z,J 



(17) 



Figure 3: The transfer function T{k) vs. y = k r/,„ for Fermions and Bosons 
decoupling ultrarelativistically and for particles decoupling non-relativistically 
(Maxwell-Boltzmann statistics). T(y) shows little variation with the statistics of 
the DM particles. We see that T(k) decays for increasing k with a characteristic 
scale ~ l/r/,„ ~ kfs which is the free-streaming wavenumber [see eq. )18H . 



Here the subindex eq refers to equilibration, the beginning of 
the MD era, I + Zeq - 3200 and T{k) is the transfer function 
which takes into account the evolution of the density fluctua- 
tions during the matter dominated era. T(k) enjoys the proper- 
ties r(0) = 1 and T{k — > oo) = 0. Namely, the transfer function 
T{k) suppresses the large k (small scale) modes. 

It is convenient to introduce the dimensionless variable 



y = k run where r,,„ = — ^ - — 



(18) 



^ V2 

V3 kf, 

Ifs and kfs stand for the free-streaming length and free- 



streaming wavenumber respectively, 

riin = 2 -^l +Zeq 



3 Ml, 



(19) 



Mo y/^DM Qp) 

Hq stands for the Hubble constant today and Mpi for the Planck 



mass. 



Ho = 1.503 10""" eV , Mp, = 2.43534 10"* GeV . (20) 

r/,„ is the characteristic length scale in the linear regime. 

We plot in fig. [3] the transfer function T{y) for Fermions 
(FD) and Bosons (BE) decoupling ultrarelativistically, and for 
particles decoupling non-relativistically [Maxwell-Boltzmann 
statistics, (MB)]. We see from fig. |3]that the transfer function 
Tij) decreases by an amount of order one for y increasing by 
unit. Therefore, T{k) decreases by an amount of order one when 
k increases by an amount of the order of the wavenumber kfs 
[see eq.([r8lll. As we see from fig. |3] T{y) shows little variation 
with the statistics of the DM particles. 



4.1. Obtaining the primordial phase density by matching the 
observed and the theoretical surface density 
We match in this section the observed surface density (Table 
1) with the theoretical surface density computed here below. 
This gives as a result eq.(l36b which determines the primordial 
phase density. 

We first compute theoretically the linear profile and the halo 
radius obtained from it to obtain the theoretical surface density. 

The linearized density profile follows from the Fourier trans- 
form of the density fluctuations today 111 j 

Piinir) = T-V r kdk sin(^r) A{k,z = 0) , (21) 
27!-- r Jo 

More explicitly, from eq.dTTb the linear profile density pii„{r) 
turns to be the Fourier transform of the density fluctuations 
MI<:,Zeq) by the end of the RD era times the transfer function 
T{k) with the result: 



Piin(r) = 



108 V2 ^DM M-pi 



xbo bi 



5n 

7,2-«,/2 

fo 

n,/2 

r, r 

lin 



//() 



(l+z,„)A|Ao| 



I dy N(y) sinly — | 
Jo \ run j 



(22) 



where |Ao| stands for the primordial power amplitude, is the 
primordial spectral index, k^ is the pivot wavenumber used by 
WMAP to fit the primordial power, k^q the horizon wavenumber 
by equilibration and 



N{y) = y 



nJ2-\ 



In 



Co y 

^eq rUr. 



T(y) 



(23) 



The numerical values of the cosmological parameters entering 
in eq.(l22li are 



|Ao 



: 4.94 10" 



Gpc" 



n, ^ 0.964 , ko^2 Gpc 



CO ^ 0.11604 



(24) 



All fluctuations with k > k^q that were inside the horizon by 



equilibration are relevant here ll 111 . This introduces in eq. (l22l) 
the comoving horizon volume by equilibration 



3 



b\ bp 



(25) 



where bi ~ 1 [actuaUy, Z?i = 1 in UM\ and bp ^ 3.669 lO^^ 
ill- 

The initial power fluctuations are multiplied by a Gaussian 
random field g(k) with unit variance 



<g{k) g\k') >^6{k-k') 



(26) 



which describes the random quantum character of the primor- 
dial fluctuations. 

Each realization of the random field g{k) with unit vari- 
ance and zero average produces a DM configuration in the 
linear regime (a 'galaxy'). The simplest one is obtained for 
g(k) = 1. The presence of g{k) will produce a large variety 
of non-spherically symmetric galaxy configurations in a large 
range of masses and sizes. For simplicity we restrict ourselves 
here to the case g{k) = 1 and leave the inclusion of g(k) + 1 
to future work. The profile punir) [with g(k) - 1] bears the 
universal properties of the galaxies, that is to say, the general 
properties common to all (or most) galaxies. This is why the 
simplest profile pii„{r) [with g{k) = 1] is very appropriate and 
useful to extract these universal properties. 

From the results eqs.(l22li-(l24li we compute and analyze the 
surface density and the density profile. 

We see from eq.(l22li that p;,„(r) decreases with r having r;,„ as 
characteristic scale since it depends on r/r/„, being the Fourier 
transform of a function of y that decreases with unit character- 
istic scale in y [see fig. |3|. 

We plot in fig. IB.BI the ratio 

Piinir) ^. . ir^(y) sin(yy) dj r 

— 77^ = 'i'iy) = 7S where y= — , (27) 

PiinW y yN(y)dy nb, 

for Fermions (FD) and Bosons (BE) decoupling ultrarel- 
ativistically and for particles decoupling non-relativistically 
[Maxwell-Boltzmann statistics (MB)], ^(y) mainly depends on 
known cosmological parameters and fundamental constants and 
has a weak logarithmmic dependence on the DM particle mass. 

We compute theoretically the surface density from the linear 
profile eq.(l22li and the halo radius eqs.(fT9]l and (IB. Il l obtained 
from first principles evolving the density fluctuations since the 
end of inflation till today. For 



f^Qlin = ro p//„(0) 



(28) 



log'io^ Vs. logioOT,; 




Figure 4: The common logarithm of the primordial phase-space density q = 
(Z Qhalo)/0^s'V)'* vs. the common logarithm of the virial mass of the galaxy 
nty = M,.,„n;/[10" Mq]. q follows by solving the trascendental equation (36) 
obtained from first piinciples for the observed values of i^oobs in Table 1 . 



we find from eqs. (fT9]l-(l22l) and (IB. lb . 



1 08 V2 

5;r 



kjMj, 
Ho a 



xbo bi 



i'J6 



J N{y) dy . (29) 



The DM linear profile eq.(l22ll decreases with the character- 
istic length rii„ which should of the same order of magnitude 
than the halo radius rp in the empiric density profile eq.([T]l. 
We define the coefficient a as a s rii„/rp and determine it 
by fitting the linear profile to the Burkert profile in Appendix 
Appendix B The value of a turns to be between 0.4 and 0.8 



depending on the DM particle statistics (see Table B.l). 

Using the numerical values of the parameters eqs. ( l20b and 
( l24l i. this theoretical formula takes the form 



tioiin = 391.064 



(MeWy bi 



where 



qp = 



Qp 



3 „8 



(keV)4 



y N(y) dy , (30) 



(31) 



and 



A?(y) = /^/2-i \n(dpql r) T{y) , do 



556.6976 .(32) 



From now on we use the dimensionless primordial density qp. 

We identify the observed surface density fioohs with the theo- 
retical value obtained in the linear approximation //q ;,«. We thus 
obtain the foUowing trascendental equation in the variable qp: 



qp 



f 

Jo 



y N(y) dy ■ 



bi 391.064 (MeVy 



(33) 
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We compute the integrals in eq.(l33]l using N{y) eq.dSTTi [i. e. 
the transfer function T{-y)] from the solution of the linearized 
Boltzmann-Vlasov equation obtained in | IJ, ^: 



y"''^ T{y) lnydy=l.3U5... , 



/-*oo 

I y"''^ T{y) dy = 2.6655 . . . 
Jo 



(34) 



and hence. 



y Niy) dy = 18.1661 (l + 0.04891 In^,,) . (35) 

These values correspond to fermions decoupling ultrarelativis- 
tically at thermal equilibrium. Bosons and particles obeying the 
Maxwell-Boltzmann statistics yield similar results as one sees 
fromfigs.[3]and lBl3] 

For fermions decoupling ultrarelativistically at thermal equi- 
librium, eq.(l33]l takes then the form: 



(1 +0.0489106 Inqp) 



1 



bi 10326 (Me V)3 



(36) 



where we used the numerical values in eqs.(l24li and ( l35T l. The 
value of ~ 1 which provides the best fit to the halo radius is 



b\ ^ 0.8 (see appendix Appendix B i. 

We proceed now to solve numerically eq.(l36]l to obtain the 
primordial phase-space density qp for different values of fioots 
given in Table 1 . 
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Figure 5: The common logarithm of the self-gravity decreasing factor Z com- 
puted from eq. B7t with qp solution of eg. 061 [fig.|4] vs. the virial mass of the 
galaxy m,. = M,,;,-,v,//[10"Mo]. 



5. The DM particle mass and the decoupUng temperature 
from the galaxy surface density 

We plot in fig. |4]the solution of eq.(l36]l, qp vs. m,,. From 
eqs.(|6]l and ( ISTT i qp can be expressed as 



ZQ. 



halo ^3 8 



(keV)4 



(37) 



Therefore, taking the observed values of the phase-space den- 
sity Qhaio (fig- ID yields the factor Z as a function of the virial 
mass m,, [eq. ©J. 

In fig. |5] we plot logjg Z and logjQ 2/„'/o vs. m,,. We see that 
Qhaio decreases with m,, while Z increases with m,, in such a 
way that the product Z Qhah is roughly constant. Moreover, 
as follows from eqs.© and (fTlT l Z Qhaio gives the DM particle 
mass 



m in keV vs. 



1.7 



nT = 49.032 Z Q, 



halo 



(38) 



We notice in fig. |5]that the factor Z changes by about two orders 
of magnitude 

2.9 10^ < Z < 5.4 10^ , 

over a large range of values of the virial mass. The variation of 
Z is relevant in the context of galaxy formation but not for the 
particle DM determination. 
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Figure 6: The DM particle mass m in keV vs. the virial mass of the galaxy 
m,, = M,.„-,„;/[10"Mo] following eq. )391 with the values of cjp solution of 
eq. )36) . We see that the DM mass in exhibits the same variation with m,, than 
the surface density /Joobs in fig- ffl The precision in the observations of the 
surface density /io translates on the precision of the DM mass m. 
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We obtain the DM particle mass m from eqs. (fT2l i-(fT3]l in 
terms of the invariant phase-space density Qp-. 



m — mo 



Ql 
keV 



mo qp 



mo = 



6. Non-universal structural galaxy properties. 

The characteristic length of the linear profile r/,„ eq.(fT9ll takes 
the following form in terms of qp eq. On i: 



keV , (39) 



where 



mo = 2.6462 keV/c^ 
mo = 2.6934 keV/c^ 



for Dirac fermions 
for scalar Bosons. 



(40) 



The numerical coefficients here correspond to ultrarelativistic 
decoupling at thermal equilibrium. For decoupling out of ther- 
mal equilibrium the coefficients are of the same order of mag- 
nitude ifioli . 

In fig. |6]we plot m according to eq.(|39]l with the values of qp 
solution of ea.(l36]l (fig. |4]l and fioobs given in Table 1. 



0.8. More 



We find m about 2 keV (up to ±10%) for bi 
generally, m is in the keV scale for bi ~ I. 



The variation of the observed surface density /loohs with the 
core radius ro (fig. [U is similar to: 

• (a) the variation of the DM particle mass m displayed in 

fig.lSl 

• (b) the variation of the primordial phase-space density qp 
in fig. m 

• (c) the variation of the density contrast in fig. |9l 

Therefore, the precision in the observations of the surface den- 
sity yUo translates on the precision in the evaluation of the DM 
mass m. 



From the solution for qp eq.(l36]l and fig. |4]we can also com- 
pute the number of ultrarelativistic degrees of freedom at de- 
coupling gd and therefore the decoupling temperature Tj which 
is a further relevant characteristic magnitude of the DM particle. 
For Dirac fermions decoupling ultrarelativistic ally at thermal 
equilibrium the number of ultrarelativistic degrees of freedom 
at decoupling can be expressed from eq.lflSll as 

gd^ 1365.5 qj. (41) 
And from fig. |4l 

0.14 < qp < 0.3 , 0.61 < qj, < 0.74 . (42) 

We thus find that gd is in the interval 

833 < gd < 1010 , 

which correspond to decoupling temperatures Td [eq.(fT4l)l 
above 100 GeV. 



r,i„ ^2\.\qp' kpc 



(43) 



In fig. |7]we plot r/,„ from eq.(l43ll and a ro from the data in 
Table 1 as functions of m,,. 

The halo radius in the hnear approximation is given by r/,>, = 
a ro which for DM Dirac fermions becomes 



ro = 



rib 



0.688 



= 30.7 qp ' kpc , 



(44) 



where we used a - 0.688 obtained in appendix Appendix B by 
fitting the Burkert and linear profiles. 

Using the range of values of qp eq. (l42l i obtained by solving 
eq.(l36ll yields 

46 kpc < ro < 59 kpc . 

which is in the upper range of the observed ro values in Table 
1. Namely, the linear approximation for the halo radius give 
values above or in the range of the observations. 

The total mass of the galaxy Mgai follows by integrating the 
density profile eq.(l22li. We find 



Mg„, ^ 20 rl pii„(0) = 20 r^ jjoohs 



(45) 



In fig. [8] we plot MgailMyiriai vs. m,, where the observed m,, and 
Myiiiai are defined by eq.(|2]i. 

We see that the ratio MgailMyirUii turns to be in the interval. 



0.12 < 



gal 



< 5. 



The contrast density, that is, the ratio between the maximum 
DM mass density p/m(0) and the average DM mass density poM 
in the universe results 



contrast s 



PlnM 
pDM 



with Pom - ^dm Pc and Q^m and pc given by eq.(fT6]l. pimiO) 
is given by eq.(l28Tl as 



PliniO) = 



P-0 lin 

ro 



We plot in fig. |9]the contrast density 

/^O lin 



contrast 



(46) 



^DM Pc ro 

As seen from fig. |9] the ratio obtained is between 1/3 and 1/2 
of the observed value ~ 3 x 10^ in 13311 . The values obtained 
are below the observed values because the linear halo radius 
'"0 - rii„/Q.6di^ is larger than the observed halo radius ro and 
the density contrast goes as 1/ro eq.(l46]l. This property shows 
again that the larger and more dilute is the galaxy the better is 
the linear approximation for non-universal quantities (see Table 
3). 
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Figure 7: The halo radius in the linear approximation ro in kpc from eg. 1441 in 
broken green line, the halo radius ro in kpc from the data in Table 1 in solid red 
line vs. the virial mass of the galaxy m,. = M,,„ ,n//[10" Mq]. ("o computed from 
first principles approaches asymptotically ro for large galaxies. 



iogiolMgai/M.irial] VS. OT,t 
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Figure 8: The common logarithm of the predicted total mass of the galaxy Mg^i 
given by ea. )45t . divided by the observed virial mass M,,iriai vs. the observed 
virial mass of the galaxy Hi,, = M,vr;a//[10"Mo]. The ratio M^dlMyiriai turns 
to be in the interval 0.12 < Mgai/My/riai < 5.0. Notice that the discrepancy of 
Mgai with Myiriai is irrelevant to the determination of the DM particle mass. 
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Figure 9: The ratio piin{0)/pDM between the maximum DM mass density p/,„(0) 
and the average mass density in the universe poM vs. the virial mass of the 
galaxy my = M,.,>,o//[10"Mo]. The ratio pHn(0)lpDM turns to be between 1/3 
and 1/2 of the observed value ~ 3 X lO' (33ll . 



Notice that we consider the whole range of galaxy virial 
masses going from 5 to 300 xlO" M^. Universal quantities as 
the surface density stay constant up to +20% within this wide 
range of galaxy masses. 



It is relevant to evaluate the halo velocity given by eq.( IA.4l ) 

v^haio = 2.316 G flora ■ (47) 
Using eq.(l44l) this equation becomes 



v'^haio = 6.705 (km/sec)^ q ' . 



(48) 



From Table 1 the input observed surface density takes the value 
juo ^ 6000 MeV^ . (49) 
Eq.(l48Tl thus becomes 



h 201 
Iv^haioiin = — r km/sec 



(50) 



The obtained range of values of qp eq.(l42li yields qp ^ 0.77 and 



halo I in 



260 km/ sec 



(51) 



This value is to be compared with the values arising from //q 
and eq.(l47li and the observed values ro in Table 1. 



(52) 



79.3 km/sec < Vv^haio < 261 km/sec . 
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The halo central density in the linear approximation is given 
from eqs. i44\ and ( |49] l by 



pou„ = ^ = 2.90 10 

'"0 



-25 



1p 3 

cm-^ 



Using the range of values of qp eq. (l42l i obtained by solving 
eQ.(l36ll yields 



1.18 10 



-25 



<po < 1-52 10"^^ , 
cm-* cm-* 



which must be compared with the observed values in Table 1 . 

We see that the linear approximation produces halo central 
densities smaller or in the range of the observations and halo 
velocities larger than the observed ones by a factor of order 
one. 

Clusters of galaxies exhibit halo radius ro about 210 kpc Isfc 
well beyond the linear halo radius ~ 50 kpc. Hence, clusters of 
galaxies cannot be described by the initial conditions used here. 
We restricted here for simplicity to initial primordial conditions 
with g(k) = 1 . Chosing general random fields g{k) + 1 fulfilling 
eq.(l26ll will provide general configurations with a large range 
of masses and sizes. Each realization of the random field gijbj 
produces a possible galaxy configuration. The factor g(k) mul- 
tiplies the transfer function and therefore is to be added in 
the r. h. s. of eqs.(fT7li. (l23T l and ( l32l i and inside the fc-integrands 
[r. h. s. of eqs.dH, (|27]i, dig, (O, (O, and We 
leave the inclusion of the factor g(k) to future work. 

In summary, the solution of the linearized Boltzmann-Vlasov 
equation presented here provides a satisfactory picture of the 
general galaxy properties. Although nonlinear effects and 
baryons are not taken into account, the linear description pre- 
sented here qualitatively reproduces the main non-universal 
and general characteristics of a galaxy summarized in Table 2. 
Moreover, the agreement is even quantitative (approximatively) 
for the Unear halo radius ro, the galaxy mass Mg,,/, the linear 

halo central density po and the halo velocity v^imioim com- 
pared to the respectived observed values in the limiting case of 
large galaxies (both ro and Mgai large). The agreement is very 
good for universal galaxy quantities as the surface density and 
the density profile as discussed above. 

7. The linear profile density: cores vs. cusps 

The properties of the linear profile density p/m(r) depend cru- 
cially on the free streaming length r/„, and therefore on the mass 
of the DM particle as we discuss here below. 

We find from eqs.(l28Tl, ( l30l l and (l35T l for the linear profile at 
the origin 



for DM particles decoupling ultrarelativistically at thermal 
equilibrium with mo given by eq.(l40t. Then eq.dSSll can be writ- 
ten as 



p,,.„(0) = 0.472977 10 



-25 /_^\ 



1.976 



\keV 



— 

keV/ 



cm^ 



(55) 



where we used the numerical values from eqs. (l24l) and ( l40b 
and the conversion of units: 



(MeV)3 
kpc 



= 0.1483698 10 



-26 



cm' 



For the DM particle mass value m ~ 2 keV found in the pre- 
vious section, p/,„(0) from eq.dSSll is two to three times smaller 
than the observed values (as it is the contrast density, discussed 
in the previous section). This is not surprising since p/,„(0) is 
not an universal quantity and given the linear approximation of 
our theoretical computation. 

We derive in Appendix Appendix C eq.( |C.9t the density 
profile behaviour for r > r/,„ where r/,„ is given by eq.(l54]l: 



r>il,„ 

Plinir) 



X 



10 



-26 



g / 36.447 kpc 



cm-' 



1.482 



In 



7.932 Mpc^ 



1 



m 
keV/ 



(56) 



It should be noticed that this behaviour has only a mild loga- 
rithmmic dependence on the DM particle mass m. The scales 
in eq.(l56]l only depend on known cosmological parameters and 
not on m. 

We present in this paper clear evidences for a DM particle 
mass around 1 keV. In addition, one can wonder what is the 
shape of the linear profile and the value of the linear density at 
the origin for a typical hundred GeV wimp. 

Since wimps are supposed to decouple non-relativistically, 
eq.(fT2li does not apply to them. For DM particles decoupling 
non-relativistically Qp is given by ifioll 



Qp 



^DM Pc 



gd (m Td)^ nonrelativistic decoupUng .(57) 



For a 100 GeV wimp decoupling at the typical temperature 
Tdwimp - 5 GeV, we find from eqs. (13111 and ( l57b 



wimp 



0.3166 10 



21 



(58) 



where we used that gd - 80 at such decoupling temperature 
1 2511 . We then find from eq.(l53]l the central density value pii„iQ) 
for such value of qp: 



PiiniO) = ^ = 336.688 b, q}' 
We use from eqs.(|39]) and ( l44b that 



([1 +0.0489106 Inqp] 



(MeV)3 



Pto,(0),w,«p- 1.773 x 10-'^ ^ 
,(53) 



cm-' 



(59) 



m 

1p = \ — 

m 



riin - 77.23 kpc 



(54) 



kpc This value for the wimps density profile at the origin turns to 
be larger than the observed values by eleven orders of magni- 
tude. This result gives a further clear evidence in favour of DM 
particles at the keV scale which reproduce very well both the 
surface density and the density profile at the origin. 
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Observed Values 


Linear Theory 


Wimps in linear theory 




5 to 52 kpc 


46 to 59 kpc 


0.0045 pc 


Po 


1.57 to 19.3 X 10-^^ A 


1.49 to 1.91 X 10-^^ A 

cm-' 


1.773 X IQ-^* A 

cm' 




79.3 to 261 km/sec 


260 km/sec 


0.0768 km/sec 



Table 2: Non-universal galaxy quantities from the observations (Table 1 combined with the virial) and from the linear theory 
results. The coiTesponding dark matter particle mass is plotted in fig. |6]and is in the range 1.6-2 keV. The larger and less 
denser are the galaxies, the better are the results from the linear theory for non-universal quantities. The last column corre- 
sponds to 100 GeV mass wimps. The wimps values strongly disagree by several orders of magnitude with the observations. 



The free-streaming length r;,„ is the characteristic scale 
where pz,„(r) varies (see fig. IB. 13b . This length is of the or- 
der of hundred kpc for keV mass scale DM particles as shown 
by eq.(l54]i. For a hundred GeV wimp decoupling at TdKimp - 5 
GeV we find from eqs.(l43Tl and dSSl ) 



runim. 



Wimp 



100 GeV, Td„i,„p = 5 GeV) = 0.0031 pc = 639 AU .(60) 



Therefore, with such small r/,„ for wimps we can use for all rele- 
vant galactic scales the asymptotic behaviour of pu„{r) eq. dC.SI l 
valid for r » r/„,. That is. 



Plinir) 



wimp 



r>0.003pc ,A 

~ = 0.8064 10"'"* 



cm' 



0.0031 pc 



1 -1-0.04616 In 



0.0031 pc 



61) 



This profile clearly exhibits a cusp behaviour for scales Ipc > 
r > 0.003 pc. Notice that this asymptotic formula eq.(|6TJ ap- 
proximatively matches around r ~ 0.003 pc the value of the 
wimp profile at the origin eq.(l59]l. 

In summary, the linear profile p/,„(r) eq. (l22l i exhibits a cusp 
around the origin for a wimp DM particle and a core behaviour 
at r = for a keV scale DM particle mass. 

To illustrate the use of the linear approximation we plot in 
fig. [10] the density profile p;,„(r) according to eqs.(l22b and 
(l27b for DM particle masses of 1 and 2 keV and the Burk- 
ert density profile for the largest galaxy ro = 51.8 kpc and 
p(0) = 1.57 X 10-^-*^ ^ in Table 1. The linear- profile best 
follows the Burkert profile for a DM particle mass slightly be- 
low 2 keV in agreement with fig. |6] 

We display in fig. [TT]the linear density profile for 100 GeV 
wimps and the NFW profile for the largest galaxy in Table 1. 
The linear density profile for 1 -2 keV particles in fig. [TO]and the 
linear density profile for wimps in fig. [TT] practically coincide 
for r >30 kpc while they strongly differ at smaller scales (r < 
30 kpc). The keV mass linear profile exhibits a core like the 
Burkert profile while the wimp linear profile exhibits a cusp 
like the NFW profile. 

We can also evaluate the halo velocity for wimps from the 
general formula eq.lfSOll and the value of qp Ki„,p eq.dSSll. We 
obtain 



v^haio lin wimp = 0.0768 km/sec 

three orders of magnitude below the observed halo velocities 
eq.(l52li. Recall that keV scale DM particles yield a halo veloc- 
ity eq.dSTI) of the same order of magnitude than the observed 



-24.5 



Burkert profile 
Plinir) for m = 1 keV. 
pii„{r) for m = 2 keV 
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Figure 10: The common logarithm of the linear density profile p/,„('") according 



to eqs. 122) and 127) in g/cm^ 



kpc and the Burkert profile eq.(T). The 



Burkert profile is plotted with red crosses for the largest galaxy in Table 1 with 
ro = 51.8 kpc and p(0) = 1.57 X 10"^^ Notice that the agreement of the 
linear profile with the Burkert profile is best for a DM mass slightly below 2 
keV. 
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halo velocities. Therefore, keV DM particles may solve the 
problem in the halo velocities recently noticed by 12711 for the 
bullet cluster when CDM wimps are used. 



-22.5 



' fJFWprofi'le 
logiop/,„(r),„„,,, vs. rinkpc 



8. Conclusions 

Dark matter is characterized by two basic quantities: the DM 
particle mass m and the number of ultrarelativistic degrees of 
freedom at decoupling gd (or, alternatively the decoupling tem- 
perature Td). By computing the evolution of the density fluctu- 
ations since the end of inflation till today we obtain the density 
profiles and theoretical relations between m and gd involving 
the observable densities poM and fio eqs.(fT2b. ( fTsT i and ( |29] l. In- 
serting the observed values of poM and jUq in these theoretical 
relations yields m, gd and Qp eqs. (I39ll-(l40t and i4T[ . respec- 
tively. 

We derive the values of the halo radius ro and the values of 
the factor Z characterizing the reduction of the phase-space den- 
sity since the equilibration time till today. For these results we 
use the observed values of the halo phase-space density Qaa/o- 

From the observed values of the surface density we present 
here clear evidence that the mass of the DM particle is about one 
or two ke V. Evidence based on the phase space density pointing 
towards a DM particle mass in the keV scale was presented in 
refs. 

In addition, one can wonder what would be the results for 
heavy wimps. For example, for wimps at = 100 GeV the 

characteristic scale r;,„ eq.(fT9]l takes the value given by eq.(l60ll. 
For such small r/„, the linear profile piinif')wimp appears as a 
cusped profile when observed at scales from 0.003 pc 1 pc as 
shown in fig. [TT] Cusped profiles are thus clearly associated 
to heavy DM particles with a huge mass m^imp well above the 
physical keV scale while cored profiles are associated to DM 
particles with mass in the keV scale. 

Notice that the linear density profile turns out to be cored or 
cuspy depending on the DM particle mass m. For m ~ keV 
the resulting linear density profile is cored as depicted in fig. 
[TOl while for m > GeV the linear density profile turns to be 
cusped as shown in fig. [TT] Figs. [TOlfTTI show that the linear 
density profiles for a 1-2 keV DM particle are similar to Burkert 
(within a factor 2-3, irrelevant for the aims of this paper) while 
for a wimp DM particle, the linear density profile is similar to 
a NFW profile. 

Interestingly enough, it is possible to derive the value of the 
surface density fio from CDM simulations. Values of the prod- 
uct ps from NFW fits to CDM simulations for galaxies were 
reported in 12211 . From these values of Tj we can derive the 
surface density jUq, since yuo = po '"o - 25 p^ with the result 



10^ Mo/pc^ 



(62) 



[Notice that p, in ll22ll dififers by a factor four from eq.(fT0b1. 

We see that the surface density from CDM simulations is five 
orders of magnitude larger than the observed surface density 
fioobs - 120 Mo/pc2 Snillij. 
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Figure 11: The linear density profile for 100 GeV wimps (broJcen green line) 
and the NFW profile (solid red line) for the same galaxy mass as the Burkert 
profile in fig. 1101 In all cases the densities are in g/cm^ and r in kpc. The wimps 
linear density profile follows eg. jC.8t . The wimp linear profile exhibits a cusp 
like the NFW profile. 



It is illuminating to insert in eq.(l36]l the above value of the 
CDM surface density fJ.^'^^ instead of the observed value iJQ„bs- 
This gives for the mass of the CDM particle nf'^'^ ~ 60 GeV 
which is a typical wimp mass. Therefore, the hnear approxima- 
tion also provides a consistent value for the mass of the CDM 
particles from simulations. 

These results show that our theoretical treatment captures 
many essential features of dark matter, allowing to determine 
its nature. When contrasted to the CDM surface density value 
obtained from CDM simulations (instead of the surface density 
value obtained from observations), our approach gives for the 
dark matter particle mass the typical CDM wimps mass scale 
(GeV), fully consistent with CDM simulations. 

It must be stressed that the linear framework presented here 
applies to any kind of DM particles: particles with mass in the 
keV scale reproduce all observed galaxy magnitudes within one 
order of magnitude, while wimps (m ~ 100 GeV) present dis- 
crepancies with observations of up to eleven orders of magni- 
tude. This is a robust indication that the DM particle mass is in 
the keV scale. 
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Appendix A. The average phase space-density 

Qhaio in sec. [3] follows averaging p{r) and v\^iig{r) over the 
volume. We define their average using the density p(r) eq.Q 
as weight function: 

rRrlr 



P = 



Vhalo = 



Jo'' r2 p(r) vj^i^ir) dr 



j^'" p(r) dr 



The virial radius 7?,.,, is defined by the radius where the mass 
computed from the Burkert profile eq.Q takes the value 131 1 



Ryir 

259 kpc 



(A.2) 



Therefore, knowing M{Ryii ) and /?,,,> we obtain ps and c from 
eqs. dA.SI l and ( IA.6b . For the galaxies in Table 1, we find 
23.2 kpc < r, < 62.5 kpc , 0.439 lO^^s ^/^.^3 < < 
1.087 10-25 g/cm^ , 7.64 < c < 13.1. We use the values 
of rs and ps for the larger galaxy to plot the NFW curve in fig. 
\n\ Namely, r, = 62.5 kpc andp, = 1.087 lO^^^ g/cml 
.(A.l) 

Appendix B. The Hnearized density profile. 

Both, the Burkert profile Fsir/ ro) eq.([TJ and the linear profile 
^(r/rim) eq.(l27Ti. have the same qualitative shape. To make the 
connection quantitative, we fit the linear profile with a Burkert 
profile setting 



Here, 



M{R 



Jo 



r^ p(r) dr^lnporl 



1 9 

ln(l + c) - arctanc + - ln(l + c') 



c = 



ro 



(A.3) 



EUiminating M(Ryir) between eqs.( IA.2b and ( IA.3b gives c as a 
function of po through the trascendental equation 



Po 



0.6187 10-27 jfr ln(l + c) - arctanc + ^ ln(l + c^) 

The right hand side is a monotonically increasing function of 
c. This implies that c increases when po increases. Since ro de- 
creases when Po increases (keeping constant the surface density 
/io), therefore c increases when ro decreases. For the galaxies 
in Table 1, we find 9.2 < c < 24.9 , 120 kpc < < 478 kpc, 
smaller values of c corresponding to larger galaxies. 

From eqs.dHi, Q and ((SJl evaluating the integrals in eq. dA.lb . 
we find 

p - 0.0662 Po , v^haio - 2.316 G po rl , , 
p 0.069 



Qhalo - 3^ 



For the NFW profile eq.lfTOt the virial mass takes the form 



(A.4) 



M{R,ir)^An r " 
Jo 



r p(r) dr - A n Ps r, 

Rvir 



ln(l +c)- 



1 +c 



and therefore we find for ps, 
p, = 0.30999 10-27 



cm 



ln(l + c) - 



(A.5) 



1 +c 

The observations give for c the empirical relation 113 111 

-0.13 



c = 9.7 



M{R,i,) 



(A.6) 



X - a y 



that is, r/,„ = a ro . 



(B.l) 



We look for the value of a that gives the best fit by minimizing 
the sum of squares: 

Vi'(y)-FB{ay)f forO<y<3 . 

The best fit for each DM particle statistics is obtained for the 
values of a in Table 2. We display in fig. IB.12l the Burkert 
profile Fsia y) and the Unear profiles for Fermi-Dirac, 
Bose-Einstein and Maxwell-Boltzmann statistics, respectively. 
We see from fig. IB. 121 that the profiles for Bose-Einstein and 
Fermi-Dirac statistics are better fitted by a Burkert profile than 
the profile for Maxwell-Boltzmann statistics. 

We compute the behaviour of the linear profile p/,„(r) eq.(l22ll 
for r » r/,„ in Appendix Appendix C We find that the linear 



approximation can be used for (see Appendix Appendix C i 



< r < r,„ 



where r„, 



8Mpc 



It must be noticed that the maximum radius r„„„ turns to be 
independent of the DM mass m and only depends on known 
cosmological parameters. 

We have at the origin F^(0) = -1 while *F'(0) = and 
"¥"{0) < 0. More precisely >P"(0) = -2.74 for fermionic DM 
particles. At the origin, the Burkert profile decreases with unit 
slope while the linear profile has an inverse-parabola shape. 

Galaxy profiles take an universal form when p{r)/po is ex- 
pressed as a function of r/ro. The Burkert profile is a particu- 
larly simple formula that satisfactorily reproduces the observa- 
tions. The linear profile ^{y), especially for Fermi-Dirac and 
Bose-Einstein statistics, fits very well the Burkert profile and 
therefore, ^(y) is also able to well reproduce the observations. 
Namely, the linear profile pim(r) is well appropriated for small 
and intermediate scales 

< r < r,^^ . 

This means that although the linear approximation cannot cap- 
ture the whole content of the structure formation, it can well 
reproduce universal features which are common to all types 
of galaxies as the density profile. Notice that the linear profile 
^iy) is universal as a function of y = r/r;„,. The values of r;,„ 
and piiniO) are not universal and change by orders of magnitude 
according to the halo mass. On the contrary, the surface density 
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Particle Statistics 


a 


Bose-Einstein 


0.805 


Fermi-Dirac 


0.688 


Maxwell-Boltzmann 


0.421 



Table B.3: The values of the parameter a = r;;„/ro for which the Burkert profile Fsia >') best fits the linear profile ^(y) = 
Plmir)/pim{0) , y = r/rii„. 




Bose-Einstein 
Fermi-Dirac 
Maxwell-Boltzmann 



1.5 2 2.5 3 3.5 4 4.5 5 



Figure B.12: The Burkert profile Fgia: y) and the lineal' profile *F(y) computed 
from first principles vs. y = r/r;,„ for Fermi-Dirac (FD), Bose-Einstein (BE) 
and Maxwell-Boltzmann (MB) statistics. The values of ff for each statistics are 
given in Table 2. The linear profile 'F(y), especially for Fermi-Dirac and Bose- 
Einstein statistics, fits very well the Burkert profile and as a consequence, ^(y) 
reproduces the observations as well as Fsia y). 




Figure B.13: The profiles pii„(r)/piin(0) vs. x, where .y = r/r/,„ for Fermions 
and Bosons decoupling ultrarelativistically and for particles decoupling non- 
relativistically (Maxwell-Boltzmann statistics). The bosons profile is the more 
peaked, the MB profile is the shallowest and the fermions profile is lying in- 
between. The profiles show little variation with the statistics of the DM parti- 
cles. 



defined by eq.([3]) is an universal quantity. Indeed, the the- 
oretical value of fi() that follows from the linear profile pim{r) 
eq. (l22l i can reproduce the observed values of as it has been 
shown in ll 111 . 



We use this property in section BTTI to derive the values of 
the DM particle mass m and the number of ultrarelativistically 
degrees of freedom at decoupling gj. 

As shown above the linear profile and the Burkert profile are 
the closest for r/,„ - a with a - 0.688. On the other hand, 
we know that the linear approximation always gives values for 
ro larger than the observed values, namely, the Unear approxi- 
mation improves for large galaxies ll ill . Therefore, we require 
that ru„ tends to ro = 0.688 ro for large galaxies which fixes b\ 
to be bi ^ 0.8. In any case the dependence of the results on b\ 
[which must be anyway b\ ~ 1] is quite mild. 
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Appendix C. Asymptotic behaviour of the Unear density 
profile. 

To derive the asymptotic behaviour of pii„{r) it is convenient 
to change the integration variable in eq.(l22li to 



T] = y — 

riin 

and we obtain 



y 



Plinir) 



(C.l) 



P/m(0) y'^ y N{j) dj 



Jo 

In the Umity = r/r/, 



( - I sin 7] drj 



(C.2) 



N 



we have from eq.(l32b 

¥-1 

■Co 



ln(y qp I + In 77 



where we used that T(Q) = 1 . 
Therefore eq. (IC.2b gives 

r(^^) sin(|«,) „, 



(C.3) 



where we used the formulas 113 



77'^' ' sin77t/77 = r^yj ' 

^ rj^^^ sin 77 ln77 1^77 = r^-j-j sin^— 77.5j 
/71s \ In Y 

— H COtg —775 , 

\2I 2 VJ 



(C.4) 



(/'(jic) stands for the digamma function. 

The asymptotic behaviour eq. dCJb is hence governed by the 
small k behaviour of the fluctuations A(^, Zeq) by the end of the 
radiation dominated era [see eq.(fT7Ti1. 

Using the numerical values for and co from eqs. (l24l i and 
( I32I 1 and the integral over N{y) eQ.(l35Tl, eg. dC.SI l becomes 



1 +0.168654 In 



y>i 0.411966 



y 



y 



1.482 



1 +0.0489106 In^ 



(C.5) 



We obtain for DM particles decoupling ultrarelativistically at 
thermal equilibrium using eqs.([39]l and ( l40l l. 



"^(y) = 0.770518 



77.23 kpc \'-^^^ 



\ m 



1 +0.11 1377 Inl — 



(in \ ' 
— 
keV/ 



(C.6) 



iMumerical calculation of *P(y) vs. 
Asymptotic formula for vs. y 




Figure C.14: The linear profile ^(y) vs. y computed from the numerical Fourier 
transform eq. j27t in red continuous line and computed from the asymptotic 
formula eg. jC.5t . We see that the asymptotic formula well reproduces the lineal' 
profile for y > 1 and not just for y 2> 1 . 



where we used 1 + 77j/2 = 1.482, 2 (2 + n,)/3 = 1.976. 

We plot in fig. IC.14l the asymptotic formula ea.( IC.5l l and the 
numerical Fourier transform eq.(l27b for We see that the 

asymptotic formula eq. (IC.5b correctly reproduces not only 
for y » 1 but for all y > 1 . 

We see that there exists a maximum value ymax (and therefore 
rmax) where the linear profile vanishes: 



y 



r,„„ = 7.932 Mpc . (C.7) 



where we used eqs. dC.lb . iC.Si and ( l54b . 

Notice that r„„„ turns to be independent of the DM mass 7r7 
and only depends on known cosmological parameters. 

Thus, the linear approximation can be used for 



<y < y„,ax 



<r < r„, 



where ^(y) > with y„,ax and r,„„ v given by eq.( IC.7b . 

The nonvalidity of the linear approximation beyond 8 Mpc 
reflects the fact that non-linear effects are important for small 
wavenumbers: this is consistent with the fact that we have ef- 
fectively cutted off the modes k < k^q in the linear approxima- 
tion [see eq.(l25]l and ll 111 1 as it must be. 

Combining the value of pim{0) in eqs. ( l53T l and ( fSSt with the 
asymptotic behaviour eq.( IC.5l ) yields 



Plinir) 



10 



-26 



g /42.029 kpc 



1.482 



In 



[1 +0.0489106 In^,,] . 



7.932 Mpc) 



(C.8) 
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We then find for DM particles decoupling ultrarelativistically at 
thermal equilibrium using eqs.([39]l and ( l40t . 



r>ri, 



g / 36.447 kpc 



1.482 



In 



7.932 Mpc 



(771 \ 
— 
keV/ 



(C.9) 



[34] 
[35] 

[36] 
[37] 
[38] 
[39] 



where r/,„ is given by eQ.(l54]i. It should be remarked that this 
behaviour has only a mild logarithmmic dependence on the DM 
particle mass m. The scales in eqs. (IC.8l l-( IC.9l l only depend on 
known cosmological parameters and not on m. 



As noticed in (author?) Ill 111 , the asymptotic decrease of the 
linear profile given by eq. dC.SI l is in remarquable agreement 
with the universal empirical behaviour put forward from ob- 
servations in llssll and from ACDM simulations in I37il . For 
larger scales we would expect that the contribution from small 
k modes where nonlinear effects are dominant will give the cus- 
tomary r^^ tail exhibited by the Burkert profile eq. ([TJ- 
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